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We investigate the influence of superpositional wave func- 
tion oscillations on the performance of Shor's quantum algo- 
rithm for factorization of integers. It is shown that the wave 
function oscillations can destroy the required quantum inter- 
ference. This undesirable effect can be routinely eliminated 
using a resonant pulse implementation of quantum computa- 
tion, but requires special analysis for non-resonant implemen- 
tations. 



I. Introduction 

The recent intensive interest in quantum computation 
was stimulated by Shor's discovery of a quantum algo- 
rithm for the factorization of integers [Q. Recent theo- 
retical quantum computation achievements include: uni- 
versal quantum logic gates 0, quantum error correction 
codes |M, Grover's algorithm for fast quantum searches 
H , and implementation of quantum algorithms using a 
statistical ensemble of nuclear spins at room temperature 
H- jq]. Introductory discussions on quantum computa- 
tion are presented in pUlfl]. One problem which has not 
yet been discussed is the influence of wave function os- 
cillations on the performance of quantum algorithms. To 
carry out a quantum computation, one utilizes superposi- 
tions of eigenstates of the Hamiltonian. As a result, each 
state can acquire its own phase, causing destructive and 
constructive interference. This interference is essential 
for quantum computation. 

In this paper we demonstrate the effect of wave func- 
tion oscillations on quantum interferences, using a sim- 
ple example of Shor's quantum algorithm for factoring 
the number 4. We show that fast oscillations can destroy 
required interference effects. We explain how this effect 
can be routinely eliminated for resonant pulse implemen- 
tation of quantum computation. We also point out that, 
for "non-resonant" proposals for quantum computation, 
the effects of wave function oscillations will require addi- 
tional consideration. 

II. A Model System of Four Qubits 

Consider four qubits which are prepared initially in the 
ground state. We shall use these qubits to implement 
Shor's quantum algorithm. The two left-most qubits rep- 
resent the number "x" , and the two right-most qubits 
represent "y" . Assume that initially all four qubits are 
in their ground states. The wave function of the system 



is, 



*(0-) = |0i0 o ,0i0o), 



(1) 



where 0~ denotes a moment of time, t = — e, where 
e — > 0; each subscript denotes the position of a qubit in a 
register; |0) and |1) denote, correspondingly, the ground 
state and the excited state of a qubit. 

To demonstrate how wave function oscillations can de- 
stroy the desired interference in Shor's algorithm, we use 
in this section a simplified consideration. Namely, we 
assume that Shor's algorithm can be implemented using 
three unitary transformations: The first transformation 
produces a superposition of all possible states of the x- 
register; the second transformation "computes" the pe- 
riodic function, y(x); and the third transformation per- 
forms a Discrete Fourier Transform (DFT). Next, we as- 
sume that all three transformations have infinitesimally 
short time duration, but there exist finite time intervals 
between transformations. A more realistic model will be 
discussed in section III. 

Let us now factor a smallest composite number, 4, us- 
ing a four-qubit quantum computer. (Generally, to factor 
a number TV, one should have the number of states in the 
x-register greater than A 2 but it is not important for our 
consideration.) Following Shor's method, we choose the 
only coprime number, 3. The quantum computer then 
must compute the periodic function, 

y{x) = T {mod A). (2) 

This function is, 

2/(0) = 1, J/(1) =3, 2/(2) = 1, 2/(3) = 3,... 

Then, the quantum computer must find the period of 
this function. (T = 2, in this case.) The factor of 4 can 
be found as the greatest common divisor (GCD) of two 
numbers: (z + 1) and 4 or (z — 1) and 4, where z = 3 T / 2 . 
In our case, 



GCD{z - 1,4) = GCT>(2,4) = 2, 



(3) 



gives the desired factor of 4. We shall use decimal nota- 
tion below for the states of the registers, 

|toiTOo, nino) — » \m, n), to = too + 2toi, n = hq + 1ri\. 

(4) 
In this notation the initial wave function (1) has the form: 
|0,0). 



Now let us consider the sequence of three quantum 
transformations required to compute the period, T . The 
first quantum transformation creates a superposition of 
all possible values of x, with a y-value of zero, 

tt 1 = i(|0,0> + |l,0) + |2,0) + |3,0)). (5) 

Assume that the delay time between the first and the sec- 
ond transformations is t\ ■ According to the Schrodinger 
equation, during the time interval t\ each state in (5) 
acquires a phase, —E pq Ti, where E pq is the energy of the 
state \p,q) and we put % = 1. After time t%, the wave 
function takes the form, 



^i(rn = \[\0 1 0)e- iBoon 



|l,0)e 



-4-ElQTi 



(6) 



|2,0)e 



— iE 2 0Tx 



|3,0)e 



-iE 30 Ti 



The second transformation computes the values of the 
periodic function (2). Thus, one gets from (6) the wave 
function, 



*i(T+) = ^(|0,l)e- iE °° Tl + |l,3>e 



-i.EioTi 



(7) 



|3, 1}) exp(-iE 00 ri - iE 01 T 2 )- 



(|0,3) +»|1,3> - |2,3) -t|3,3))exp(-*Biori - iE 13 T 2 )+ 



(|0, 1) - |1, 1) + |2, 1) - |3, 1)) cxp(-^^ 20 T 1 - iE 21 r 2 y 



(|0, 3} - *|1, 3) - |2, 3} + »|3, 3)) expf-iEaoTi - iE 33 T 2 )}. 

One can see the destructive impact of the wave function 
oscillations on the quantum computation. (The phase 
advance at later times t > tx+t^" does not further change 
the quantum interference.) If T\ = t 2 = 0, one gets the 
desired constructive and destructive interferences. In this 
case, the wave function (10) reduces to the following, 

tt = i(|0,l) + |2,l)+|0,3>-|2,3». (11) 

Measuring the state of the ^-register one gets two values 
for x: x\ — and x 2 — 2. The period, T, can be found 
as the ratio, D/x 2 , where D is the number of states in 
x-register fllj| . In our case, D = 4 and T = 4/2 = 2. 

If Ti,T2 =/= 0, the required quantum interference is de- 
stroyed unless the following very special conditions are 
satisfied, 

{E 20 - £oo)ti + (£21 - E Q i)t 2 = 27T7H, (12) 



|2,l)e- 



-iE 2 ori 



|3,3)e 



-iE 30 Ti 



At this point, the "natural" relation between the state 
and its phase is destroyed: the phase factor of the state 
\m, n) is not equal to —E mn r\. 

During the next time interval, t 2 , the wave function 
evolves further according to the Schrodinger equation 
and takes the form, 



(£30 - Exq)tx + {E 33 - E 13 )t 2 = 27rn 2 , 

where n\ and n 2 are arbitrary integers. For example, the 
state |1, 1), which vanishes if T\ = t 2 = 0, now survives 
as a sum of two terms, 

-|1, l)[exp(— iEfooTi - iE Q1 r 2 ) - exp(-iE 20 Ti - iE 21 r 2 )]. 

(13) 
The first term (1/2) |0, 1) in the superposition (11) is de- 
scribed now by two terms, 



*(n + r 2 ) - -{|0, 1) cxpH^oori - tE Q1 r 2 )+ (8) 



|1, 3) cxp(-i£; 10 r 1 - iE 13 t 2 )- 



|2, 1) exp(-i£; 2 ori-i£; 2 ir2)-|-|3, 3) exp(-i£;3 ri-i£:33r2)}. 

Finally, at t = T\ + t 2 one performs the DFT of the re- 
register, 



1 3 



} 2irikx/4: 



\k). 



(9) 



fc=0 



Then, the wave function takes the form, 



*(r 1 +r+) = i [(|0,l) + |l,l) + |2,l)+ (10) 



-|0,l)[exp(-iJSoo7i - iE 01 r 2 ) + exp(-i J B 20 T 1 - iE 21 r 2 )\. 

(14) 
The reason Shor's algorithm fails in this example is 
that the "like" terms corresponding to the same state, 
carry their "history" after their initialization. For exam- 
ple, the first term in (14) is generated by the transforma- 
tion sequence, 



|0,0>-|0,0)^|0,1)3|0,1), 



(15) 



while the second term in (14) is generated by the trans- 
formation sequence, 



|0,0>-|2,0}^|2,1)3|0,1). 



(16) 



One can see that the two different phase factors in (14) 
directly reflect the "history" of the corresponding terms 
(15) and (16). 



To provide a quantum algorithm for a finite time one 
must create the same phase factor at the end of the algo- 
rithm's implementation, for all like terms. The natural 
way to solve this problem is to "induce" the "natural" 
phase, —E mn t, every time the term \m, n) is generated 
in the process of quantum computation. In the next sec- 
tion we show that this can be done routinely using a 
resonant pulse implementation of quantum computation. 

III. Generating "natural" phase factors using 
resonant pulses 



Following the Lloyd's idea (Tj], we assume that all 
quantum transformations are implemented by a proper 
sequence of resonant pulses. Each pulse drives a quan- 
tum transformation between the corresponding energy 
levels in a system of weakly interacting qubits. Next, 
we assume that one keeps the stable continuous ref- 
erence oscillations, ~ exp(itt> p fct), with the frequencies, 
ujpk = E p — Ek , for each resonant transition which is used 
in the process of quantum computation. (Here E p and 
Ek are the energies of the states \p) and |fe) of the whole 
system of qubits.) The resonant pulses are supposed to 
be "cut" from the reference oscillations and to be ap- 
plied to the system of qubits with a proper phase shift, 
ip {— it < tp < tt) relative to the reference oscillations. 
This technique is routinely used in nuclear magnetic res- 
onance. 

Now we show how a natural phase is generated by 
a resonant pulse. Let's consider the state |fc) of the 
whole system of qubits. Suppose that the resonant pulse 
with frequency Lu pk drives the transition from state |fc) to 
state \p), i.e. "generates" the state \p). Neglecting non- 
resonant effects, the Hamiltonian for this process can be 
written as, 



iC p = E P C P - ine-^^Ck. 

The solutions of these equations are, 

C k (t) = C k {t )cos[rt{t - to)/2] cxp[iE k {t - t)\, (20) 

C p (t) = C k (to)sm[n(t-t )/2]exp[i{n/2-(p+E k to-Ept)}. 

At the end of the pulse: t = to + r, one has from (20), 

C fe (io+T) = \C k {t )\ cos a exp[-iE k (t + T)}, a = ttr/2, 

(21) 



C p (t - 



\C k (t )\ sinaexp[i(7r/2-y>)] exp[-iE p (t +T) 



One can see that at the end of the pulse the newborn 
state, \p), has the natural phase factor, exp(— iE p t), 
while the initial state, \k), keeps its natural factor, 
exp(— iE k t). This shows that the resonant pulse trans- 
fers the phase accumulated by the reference oscillations 
to the newborn state, \p), providing the natural phase 
factor, exp(-iEpt). The factor exp[i(7r/2 — ip)] in the 
expression for C p (to + r) is the "standard" phase shift 
which can be eliminated by choosing ip = tt/2. 

On the contrary, let us consider now the case when a 
resonant pulse is not cut from the reference oscillations, 
and has the initial phase, (po (— tt < ipo < tt), at time 
t = t{). In this case, the interaction between a qubit and 
a resonant pulse (the third term in (17)) has the form, 



n 



{exp[iu pk (t - t Q ) + iipo]\k)(p\- 



(22) 



exp[-iw p k(t - t ) - iip ]\p)(k\}. 



H = Ek\k){k\+Ep\ P ){p\ 



(17) 



In this case, the amplitude of the newborn state, \p), in 
(21) will be replaced by, 

Cp{to + r) = | C fe (t ) | sin a exp[i(7r/2-(p )]x (23) 



e l ^ kt+ ^\k)(p\ + e-* u "' t+v) \p)(k\ 



Here $7 is the Rabi frequency, and tp is the phase of the 
resonant pulse. 

Assume that at the beginning of the pulse (t = to) the 
amplitude of the state |fc), Cfe(io), has the natural phase 
factor, 



C k (t ) = \C k (t Q )\exp{-iE k to), 



(18) 



and the amplitude of the state \p), C p (to) = 0. The 
Schrodinger equation for the amplitudes C k and C p , for 
the states \k) and \p), reduces to, 

iC k = E k C k - hle^» t+ ^Cp, (19) 



exp(-i£ fc t 



iE p t) 



Instead of the natural phase (21), one has a phase which 
reflects the history of the corresponding state. Namely, 
the factor exp(— iE k to — iE p r) in (23) shows that the 
state |p) was generated from the state |fc), at t = to- 
To get the natural phase, one could connect the value of 
the initial phase, <po, to the time of the pulse applica- 
tion, to. Indeed, choosing ipo = p + ^> P kto (where ip is 
a proper phase shift: — tt < (p < tt), we obtain the ex- 
pression (21) derived for a coherent pulse. However, in 
the case of application of a non-coherent pulse, the in- 
evitable uncertainty, 5to, can cause a large uncontrolled 
phase shift, uj pk Sto- Taking into consideration that quan- 
tum computation requires application of a large number 
of pulses, one can conclude that practical implementation 



of non-coherent pulses is a complicated technical prob- 
lem. Thus, applying non-coherent resonant pulses one 
can destroy the desired quantum interference. 

The same destructive effect may occur for non-resonant 
implementations of quantum computation. Assume, for 
example, that a powerful pulse, V{t), 



V{t) 



V, if t <t<t + r; 
0, if t < t , t > to + T, 



acts on an individual qubit driving it from the ground 
state to the excited state (or, vice versa, from the ex- 
cited state to the ground state). This corresponds, for 
example, to the transition of the whole system of qubits 
from the state, \k) to the state \p). The last term in (17) 
which describes the interaction between the qubits and 
the external field, we present in the form, 



-V(t)(\k)(p\ + \p)(k\ 



(24) 



For the simplest case, V 
derive the solution, 



oo, t — ► 0, Vt — ►= a, we 



Cfc(*o") = Cfc(*o ) cosa > Cpito) = *Cfe(t ) '• 



(25) 



One can see from Eq. (25) that the newborn state, \p), 
acquires the phase accumulated by the "parent" state, 
\k) (plus the "standard" phase shift, tt/2). Again in this 
case, the phase of the newborn state, \p) (with the ac- 
curacy to the "standard" phase shift, tt/2) is equal to 
—Ekto, instead of the natural phase, —E p to- 
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Conclusions 

We have demonstrated the importance of wave function 
oscillations for finite-time implementation of Shor's quan- 
tum algorithm. These oscillations can destroy the de- 
sired interference required for quantum computation. For 
resonant implementations of quantum computation, the 
detrimental effects associated with this phase factor can 
be routinely eliminated if one can keep stable continu- 
ous reference oscillations with the resonant frequency for 
each quantum transition which is used in the process of 
quantum computation. For "non-resonant" proposals for 
quantum computation the influence of wave function os- 
cillations on performance of quantum algorithms requires 
a special consideration. 
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